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We show that the critical behavior of quantum systems undergoing a percolation transition is 
dramatically affected by their topological Berry phase 2np. For irrational p, we demonstrate that 
the low-energy excitations of diluted Josephson-junctions arrays, quantum antiferromagnets, and 
interacting bosons are spinless fermions with fractal spectrum. As a result, critical properties not 
captured by the usual Ginzburg-Landau- Wilson description of phase transitions emerge, such as 
complex critical exponents, log-periodic oscillations and dynamically broken scale-invariance. 



A fundamental aspect of the Ginzburg-Landau- Wilson 
(GLW) description of phase transitions is scale invari- 
ance, which relies on the absence of characteristic length 
and energy scales at criticality, leading to the concept of 
universality [3, For instance, near a quantum crit- 
ical point (QCP), if a physical observable O (T) trans- 
forms for an arbitrary scale transformation b > ac- 
cording to O (T) = b~ x O (b z T), then we obtain a power- 
law temperature dependence O (T) oc T x l z , with univer- 
sal critical exponent x/z. However, if scaling is valid 
only for powers of a discrete value bo, it follows that 
O (T) = T X I Z Q (InT), with Q (t) a periodic function of 
period z In fro- Fourier expansion of Q (t) yields: 
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with constant coefficients a n . Thus, the system is charac- 
terized by a family of non-universal complex exponents. 
An invariant scale bo, leading to this discrete scale in- 
variance, is found in several critical systems that either 
are out of equilibrium or have an underlying built-in hi- 
erarchical structure (for a review, see Q). 

In this Letter, we show that complex critical expo- 
nents and log-periodic behavior appear in a variety of 
disordered systems close to a percolation QCP, such as 
Josephson-junction (JJ) arrays, quantum antiferromag- 
nets (QAF) and interacting bosons. Rather than being 
related to non-equilibrium properties or to the fractality 
of the percolating cluster, in these systems the invariant 
scale bo emerges naturally in their low-energy excitation 
spectrum for certain values of their Berry phase 27tp. 

By calculating their specific heat and compressibil- 
ity at the percolation threshold, we show that, for ra- 
tional p. large clusters have the lowest excitation en- 
ergies, giving rise to usual power-law behavior below a 
crossover temperature T*, which varies in a pronounced 
non-monotonic way with respect to p (see Fig. 1). For ir- 
rational p, the low-energy properties are governed instead 
by degenerate clusters of intermediate sizes, leading to 
the breakdown of continuous scale invariance (T* — > 0). 
Remarkably, the sizes and energies of these resonating 
clusters depend solely on the continued-fraction expan- 




Figure 1: (a) In the diluted quantum system, the global phase 
(arrow) of clusters of connected grains (dark/red) coherently 
precesses due to the Berry phase 2np. (b) Strong variation of 
the crossover temperature T* below which scaling with real- 
valued exponents holds (blue points, To ~ U). 



sion of p. As a result, when p is a quadratic irrational, 
the periodicity of its continued-fraction expansion gives 
rise to an invariant scale bo and, consequently, to complex 
critical exponents. 

To introduce a model for all the systems discussed 
above [2|, consider an array of grains characterized by 
an XY order parameter = l^ol exp (i9j), with phase 
9j and fixed amplitude |Wo|. The array is diluted on a 
regular lattice of dimension d > 1, characterized by a 
quenched random-site variable ej that takes the values 
and 1 with probabilities P and (1 — P) , respectively. We 
consider the Hamiltonian 0, Q, @| : 

H = Uj2 e i ( n i ~pf~Y. ' cos & - e o)' ( 2 ) 

i ij 

where rij = —i-gjf- P can be externally controlled and 
causes the phase to precess in time according to d9j/dt = 
2Up (see Fig. 1). In JJ-arrays 0, ^ j denotes the 
superconducting order parameter, Jjj is the Josephson 
coupling, U is the charging energy and p, related to the 
AC-Josephson effect, can be changed by an external gate 
voltage, tyj can also represent planar quantum rotors, 
associated with the low-energy modes of QAF (2j, where 
p is proportional to a perpendicular external magnetic 
field. For systems of interacting bosons which can 
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be realized in optical experiments with cold atoms |8|, 
2Up corresponds to the chemical potential /i. 

The effects of percolative dilution on all these systems 
have been the subject of various experimental and nu- 
merical investigations • Here, we focus on the criti- 
cal properties at the percolation threshold P c , where the 
density of clusters with s connected occupied sites varies 
as N(s) (x s~ r , with 2 < t = d/D f + l < 2.5 and D f the 
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fractal dimension of the percolating cluster [13j- At low 



temperatures T <C \ Jij | and deep in the ordered state 
of the clean system (U < \Jij\), the relative phase be- 
tween grains inside each cluster is fixed, implying that 
the entire cluster is characterized by a global phase 
[l6| . At P c , contributions to the total specific heat of 
a single cluster arise from the coherent precession of its 
global phase, C, and from the excitations of internal col- 
lective (spin-wave) modes that change the relative phase 
between grains, C sw . As we will show below, C sw is 
sub-leading; thus, similar to the behavior in superpara- 
magnets, each size-s cluster can be treated effectively as 
a big single rotor, with the corresponding action: 
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which describes the coherent phase-precession due to 
both quantum fluctuations and the Berry phase 2irp. No- 
tice that, unlike the case of SU(2) spins, the Berry phase 
of quantum rotors has a topological character, since the 
imaginary part ^Bcrry = isp f df^ of A s is indepen- 
dent on the time evolution of 9 (r ) , enabling us to solve 
our problem using sums over winding numbers. Shifting 
the imaginary time f — > f/s in ([5]) eliminates the pre- 
factor s at the expense of a cluster size dependent tem- 
perature T — » sT and, most importantly, Berry phase 
p — » sp. This yields the free energy scaling F s (p, T) — 
s~ 1 F\ (sp,sT), from which we can derive scaling rela- 
tions for the heat capacity C s (T) = -Td 2 F s /dT 2 and 
the compressibility k s (T) = —d 2 F s /dp 2 . Here, the suffix 
1 (s) refers to quantities on a single site (single cluster). 
Thus, macroscopic quantities can be calculated by aver- 
aging over all clusters, i.e. O (p, T) = J2 S N (s) O s (p, T). 

Let us first revisit the results for p = 0, where univer- 
sal power-law behavior was previously found The 
low-temperature specific heat of a cluster is given by 
C s (p = 0) oc exp (— U/sT), i.e. the typical excitation 
energy of a cluster decreases monotonically with its size, 
e s = U/s. Then, the low-energy behavior is dominated 
by large clusters and we can replace the sum over s by 
an integral, obtaining, for T <C U, 



C(p = 0,T) cx T d/z 



(4) 



where the dynamic scaling exponent z r = Df was intro- 
duced. This result was also found in detailed computer 
simulations flTl [l8| . Consider now a finite Berry phase 
p 7^ 0. Solving the problem of a single quantum rotor 




Figure 2: Excitation energy e s as function of the cluster size 
s for the irrational p = y7. The highlighted points corre- 
spond to energetically degenerate clusters associated to the 
Diophantine approximant p/q = 37/14, and are responsible 
for a jump in the integrated density of states. The inset shows 
the spectrum for the rational p — 37/7, characterized by four 
well-defined branches. 



and using the scaling T — > sT, p — > sp derived from 
the action ©, we find the spectrum of a single cluster 
E m = U (m — sp) 2 /s, with m £ Z. Therefore, the lowest 
excitation energy is (see also [1, 0|) : 



e, = -(l-2|p.|) 
s 



(3) 



where p s = sp — [sp + ij and [^J is the integer part 
of a;, i.e. \x + is the integer closest to x, such that 
\Ps\ < 1/2. Note that p s depends on the droplet size in 
a highly non-monotonic way, reflecting the periodicity of 
the Berry phase (see Fig. 2). Now, not only large clusters 
yield small excitation energies, but also intermediate-size 
clusters with \p s \ < 1/2. 

These low-energy excitations can be described by spin- 
less fermions, yielding the well-known fermionic expres- 
sions C\ (ui, T) — (lo/T) 2 np (u>) np (—oj) and k\ (oj, T) = 
(1/T) np (cj) np (—uj) for the single site specific heat and 
compressibility, with np (w) the Fermi function. In Fig. 
3 we demonstrate this numerically by comparing C (T) 
obtained from the exact energy spectrum and from the 
fermionic expression only. We can also show this result 
analytically: Consider, for definiteness, < p < 1/2. 
The spectrum of a single rotor, including degeneracies, 
can be generated from a model of effective interacting 
fermions and bosons, with occupation numbers nt = f 
and rib = b^b, respectively 



H = U (n b + (1 - 2p) nf + pf 



(G) 



For p — 0, one recovers the N = 2 super-symmetric de- 
scription of a rotor [l{|. Expanding for p ~ 1/2, we 
obtain instead Hq ~ £i?i/ + Urib + Unl + 2einbnf. While 
the excitation energy of single bosons is U, the fermionic 
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Figure 3: Specific heat C as function of temperature T for the 
rational p = 55/89. The dashed blue line is the exact result 
and the solid red line is the fermionic approximation. The 
onset of power-law behavior is marked by T* (dotted line). 
The inset shows the log-log variation of T* with respect to 
the denominator q for a series of rationals with comparable 

vnlnpq „ = /3A2155144377 987 1 
' r \ 5 ' 13 ' 34 ' 89 ' 233 ' 610 ' 1597 I ' 



excitation energy is z\ = U (1 — 2p) <C U. Thus, at suffi- 
ciently low temperatures, bosons are diluted and the in- 
teraction terms can be neglected, implying that free spin- 
less fermions are the dominant excitations of the system. 
The limit p — > 1/2 also plays an important role in the 
excitation spectrum of droplets in the Bose-glass phase 
[fl and in the insulating phase of interacting bosons in a 
disordered chain Q. 

Let us now consider p — p/q to be rational, i.e. p and 
q are integers with no common divisors. Due to the pe- 
riodicity of the Berry phase term, the number of distinct 
values of p s is of the order of q/2, defining well separated 
branches in £ s , all decaying as s _1 , with lowest branch 
e s ~ U {qsy 1 (see inset of Fig. 2). For T < U/q 2 , the 
problem is virtually the same as for p = 0, leading to a 
heat capacity dominated by very large clusters as given 
in Eq0] with the same exponent shown in Fig. 3. 

For the compressibility, we obtain n (T) oc T d ^ z -~\ with 
a Wilson ratio (kT) /C « 0.3 for d = 2. Note, from Fig. 
3, that the crossover temperature T* changes as q~ 2 and 
is insensitive to p. Thus, systems with similar values of 
p can have very different T* (see Fig. 1). 

The natural question is: what happens in the regime 
U I q 2 <C T <C U for very large q? For irrational p, T* — > 
and this regime prevails down to the lowest energies. 
Indeed, when p is irrational, the sequence p s is uniformly 
distributed between —1/2 and 1/2 [20], i.e. there are 
finite-size droplets with arbitrarily low excitation energy 
e s . It is convenient to introduce the averaged integrated 
density of fermionic states, D (lu) = S ~ T 9 ( w — £ s); 
from the periodicity of s s with respect to p (i.e. summing 



over winding numbers) , we obtain, for 

,/saw ($X\) 



D(u/)=C(t-1)-+ J2 



(7) 
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where x± = lu/U ± 2p, C,(x) is the zeta function and 
/saw (x) is the sawtooth function, which has period 2 and 
unit jumps at every odd integer: / saw (x) = (—x + 2n) /2 
for 2n — 1 < x < In + 1. These jumps give rise to 
discontinuities in D (w) at frequencies: 
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where pj is an odd integer and qj is even. At uj <C U, the 
fractions Pj/qj that satisfy Eq. ([5]) are the ones that best 
approximate p, i.e. the Diophantine approximants, which 
are given by the convergents of the continued fraction ex- 
pansion of p [21]. Physically, the jumps are a consequence 
of the existence of a set of energetically degenerate (i.e. 
"resonating") clusters with sizes that are odd multiples 
of <Zj/2, s = (n + 1/2) qj (see Fig. 2). Summing over all 
these clusters in Eq. ([7]), we find that each jump in D (lu) 
is given by A, = qJ T ( (r) (T - 1). 

Back to Eq. [TJ we find that the regular part of the saw- 
tooth function cancels out the linear in lu term. Thus, 
the frequency dependence of D (lu) is governed by the 
successive jumps Aj at LUj of Eq. ([5]) and, consequently, 
by the sequence of convergents of the continued fraction 
expansion of p with even denominator qj. Although the 
determination of this sequence for an arbitrary irrational 
p is an outstanding problem in number theory, it is sim- 
plified in the case of quadratic irrationals, which have 
periodic continued fraction expansions [21| . Then, one 
finds that the sequence of even qj is also periodic. 

In fact, for quadratic irrationals with a single period 
a 6 Z, which are solutions of the algebraic equation 
y 2 — ay — 1 = 0, we find that if qj is even, so is qj+N, 
with N = 2 + mod (a, 2). Consequently, the distance 
between consecutive jumps is a constant in log-scale, 
In (LUj/LUj + i) = 2iVTn?/_|_, as well as the ratio between 
their amplitudes, In (Aj/A J+ i) = rN\i\y + . where y+ 
is the positive solution of the algebraic equation. Using 
these properties, we can show that D (lu) is a fractal func- 
tion with fractal dimension d w = r/2, characterized by 
a power-law decay in lu and periodic oscillations in lnw. 
Since C (T) = f dec v (lu) d (lu, T), where v = dD/diu is 
the density of states, we obtain: 



C (T) = T d/zw A (InT) . 



(9) 



where Zi r — 2Dfd/ (Df +d), i.e. zi r > z r , and A (t) is 
a periodic function of period z; r In bo = 2N In y + . For 
the compressibility, we find k(T) = T { - d / Zir ^ 1 B (InT), 
where B (t) has the same period as A (t). Thus, the sys- 
tem has complex critical exponents as in Eq. ([T]). In Fig. 
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Figure 4: (a) Frequency dependent integrated density of 
states D(u), as well as the temperature dependent (b) spe- 
cific heat C (T) and (c) compressibility k (T) for the system 
with p = \/2- Dashed lines show the underlying power-law 
behavior superimposed to the log-periodic oscillations. Inset 
is a zoom of D (to). 



scaling properties of the fermionic density of states, we 
find O sw cx OT^ for both O = C, «, with <f> = d s - 1 
(4> = d s — 1/2) for rational (irrational) p. As d s > 1 [Hj], 
it follows that <fi > 0. Since the internal modes are sub- 
leading compared to the coherent modes, the spectrum 
of a single cluster depends solely on its size and not on 
its shape, in accordance to Eq. [3l 

In summary, we solved the XY quantum-rotor problem 
at low T and close to the percolation threshold, which de- 
scribes diluted systems as diverse as JJ-arrays with d.c- 
bias voltage, canted QAF in a perpendicular magnetic 
field, and interacting bosons coupled to a particle reser- 
voir. Their topological Berry phase 2irp dramatically al- 
ters the percolation QCP, since the low-T behavior is 
governed by emergent spinless fermions with fractal spec- 
trum, giving rise to generally irregular log T-oscillations 
of thermodynamic variables. While for irrational p they 
persist to T — > 0, for rational p — p/q they occur in the 
temperature range q~ 2 < T/U < 1, which can be broad 
for <;> 1. Remarkably, for a quadratic irrational p, they 
become regular, leading to complex critical exponents. 
Our results demonstrate that the quantum criticality in 
disordered systems governed by a topological Berry phase 
is beyond the GLW paradigm of critical systems. 

We thank M. Axenovich for pointing out Ref420] to us. 
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BES, Materials Sciences and Engineering Division. 



4, we show numerical results for p = y2; we also veri- 
fied numerically that the scaling form in Eq. [9] holds for 
quadratic irrationals p with more complicated continued- 
fraction periods. For non-quadratic irrationals, our nu- 
merical calculations indicate that Eq. [5] still describes 
the critical behavior, but now A (t) oscillates irregularly 
without a well-defined period. 

The breakdown of continuous scale invariance for irra- 
tional p can be attributed to the resonating clusters with 
arbitrarily low excitation energies, as they cause jumps 
in the entire spectrum, prohibiting to replace ^ s — > J ds 
in Eq. [7] For rational p, such a replacement is al- 
lowed, leading to full scaling s — > s/b Df and to Eq. [4] 
Yet, when p is a quadratic irrational, the dynamically 
broken scale invariance is partially restored as discrete 
scale invariance. In this case, the periodic structure of 
the continued-fraction expansion of p gives rise to log- 
periodic relations between sizes s = (n + 1/2) qj and en- 
ergies ujj — 2U\p — Pj/qj\ of different sets of resonating 
clusters, establishing an invariant scale bo and complex 
critical exponents (d/zi rr ) + inn/ (N\ny + ) for C (T). 

Going back to the contribution of the spin-waves, their 
density of states is v sw oc k ds ~ 1 at P c , where the fracton 
dimension d s characterizes the spectrum of the eigenval- 
ues k 2 of the Laplacian on the percolating cluster Til, 1(| . 
From the spin- waves dispersion ft 2 k = e 2 + c 2 k 2 and the 
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